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Lew  Fre.i'uency  Instability  cf  Lisnii 
jr.r lesions  Iriven  by  Multiple  Free  Piston; 


Ir.tn;  iucticr. 


T;.e  use  -.f  i;r.plcling  lisuii  cylinders  ..r  liners  to  ccapress  plasr.a 
ind/i. r ra.pr.etic  flu;':  tc  high  energy  densities  has  recently  received  con- 


siderable attention",  .d  basic  prcble;n 


in  such  dpciar.ic  con- 


pressicn  techscicues  is  the  growth  cf  hydro i;/r.an-.ic  perturbations  at  the 
free  surfaces  of  the  liner  as  these  surfaces  accelerate  toward  the  histher 


.’.tass  density  liquid  in  response  to  tl.e  pressures  of  the  driving  fluid 

« 

(.gas  or  .magnetic  field)  and  the  payload  plasna  and  ■''or  field,  fuch 
growth  is  ccr.ple'ely  analogous  t..  the  pkencn’.encr.  cf  F.ayleigh-Taylcr  in- 


sracilitj-^  in  a gravitational  field,  and  can  result  in  both  the  disrup- 
tion cf  the  cuter  surface  cf  the  ir.plcding  liner  during  the  initial  in- 
pl^sii.r.  and  the  destruction  ^:f  ti'.s  high  ter.perature  plasma  payload  by 
inner  surface  material  near  peal;  cempression.  The  latter  effect  can 
elimi.nate  the  liner  implcsicn  technique  as  a method  cf  achieving  fusion 
plasr.as;  v/hile  the  form.er  process  severely  lim.its  the  possibility  cf  con- 
trolling and  recovering  the  liner  ki.netic  energy  which  is  im.pcrtar.t  in 
fusion  reac‘‘or  applications". 


*l.at  the  Payleigh-Taylcr  instability  at  the  inner  surface  near  pea.-: 
Note:  Manuscript  submitted  May  31,  19T7 


1 


’un  i:e  o Ii;:;inatoi  L-y  jpirmiag  the  li;:er.  "faring  ii;.plL.:icr:, 
the  aein.uthal  Jcee.i  i-a"  the  inner  surface  i:;ateriai  ittcreasej , iue  ti. 
eerxervatii-;'.  ei"  ainguiar  i.-xr.entui::,  euch  tint  the  een.tripetal  acceleraticn 
. vercor-eo  the  raaial  acceleraticn  (r)  ti.  reverse  the  effective 
gravity  at  the  interface  (r  - v^/r)  in  favcr  i,f  stability,  cuch  rotation, 
hcv.-ever,  doec  net  rer.cve  the  Rayleigh-Tayicr  instability  at  the  cuter  sur- 
face v/hen  the  liner  is  accelerated  initially  by  a gas  .;r  magnetic  pressure 


and  similarly  when  i' 


ilerated  after  rebcuniinr  frc:r.  pea^i  pap'lot 


cci.'.pressi^n . 

To  eliminate  the  Kayleigi.-Taylcr  instability  at  the  free  ^uter  sur- 
face, it  •.;as  suggested'  tint  the  free  cuter  surface  itself  be  eliminated. 
The  liner  •.•rculd  be  driven  by  free-pist^ns  in  ci.ntinuous  contact  v;ith  the 
fluid,  '..•ith  the  pistons  providing  a stiff  interface  between  the  high 
uensity  lipuid  liner  and  a lev;  density  fluid  drive,  such  as  high  pressure 
gas.  The  criminal  concept  v;as  t,.  use  a fiat  .iisc  or  annulus  displaced 
axially  to  dcive  the  fluid.  The  first  e:-:perir;ental  arrangem.ent , however, 
consisted  of  a plurality  of  free-pistor.s  displaced  radially  inv;ard  to  in- 
ject fluid  int>-  a central  implosion  cham.ber.  This  is  shewn  schematically 
in  Fig.  1.  Experim.ental  tests  with  such  systems  dem.onstrated  that  the 
cuter  surface  instability  was  indeed  subdued,  with  repetitive  im.plosicn- 
reexpansion  cycles  observed.  It  was  pointed  out“,  however,  that  the 
asim.uthal  array  of  free-pistcr.s  should  be  subject  to  Rayleigh-Taylcr  in- 
stability for  ’wavelengths  larger  tl'.an  the  pisten  diameter.  E.x pe rimer. - 
‘ally,  asi.muthal  asymmetries  in  the  radial  pi.sitions  of  the  pistons  are 
bserved.  '.Vhile  'various  pr..cesses  car.  contribute  to  such  asimemetries , 
it  is  useful  to  consider  first  the  grtw'hi  of  p>  siti^  nal  variati- ns  due 
hydro.i,v:-.3mic  instability. 


o 


TI.o  aiia2y.:ij  . i'  the  jcctii^n  ij  Jir.llar  ii;  jpiria  U the 

Rayleiph-Taylcr  aualyjie  but  differs  in  that  the  interiace  be- 
•v;ee;.  the  heavy  and  liitht  fluids  cuneitte  a Jtiff.  thin  plate  .iio- 
p laced  mrv.'.al  to  itc  curi’ace.  A auececcion  of  oucii  plateo  oi’  piotono  ic 
dioplaced  i;.  a periodic  manner  to  pia^vide  the  in.itial  condition  of  the 
o..-undary  bcf.veen  the  tv/c  fluidc.  In  the  clacoical  Rayleigh-Taylor 
arialyoio,  f'.r  each  tr.i.  de  tiie  interface  io  a cine  cr  cosine  v/ave.  In  the 
preser.t  v;crl. , the  interface  at  any  tir.e  appears  as  a train  ^ f square 


■ ,i.-c..ei 

From  the  ab>^ve  discussion,  it  is  clear  titat  the  fundan.ental  prcbler. 
is  that  of  the  hydro d.vr.amic  or  Rayleigh-Taylor  instability  ^-f  a surface 
comprised  ^ f finite  elements.  To  isolate  the  main  physics  of  this  problem 
in  its  most  simplified  f rn.,  v;e  neglect  curvature  and  consider  the  situ- 
ation depletes  in  Fir.  2.  That  is.  v;e  treat  an  i:-.finite  fluid  slab  of 
iepth  d in  a constant  vertical  travitational  field  g.  The  upper  surface, 
s - i.  is  taken  to  be  a rigid  ■.■.•all.  The  lotver  surface  s = -d  is  taken 
to  be  comp,  sed  of  fir.ite  rigid  ele.ments  of  length  £ in  tl.e  x direction. 
This  then  represents  a variation  of  the  classical  Rayleigh-Taylor  in- 
stability problem  v;here  for  a continuous  free  surface  at  z = -d  the  sur- 


face perturbation  growth  rates,  uj,  satisfy^ 


- k.c  tanh  (ko 


where  k is  the  horizontal  wave  number  of  ti;e  disturb.ance  k ,=  2t:  'X  and  X 
is  its  •.^fave length.  The  sho'rter  wavelength,  dis'urbances  grow  faster  than 


the  longer  ones. 
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For  Che  finite  element  surface  configuration  of  Fig.  2,  Che  situ- 
ation is  drastically  changed.  Now  a constraint  is  imposed  on  the 
Fourier  harmonics  of  the  surface  displacement  to  insure  that  the  re- 
sultant surface  is  always  consistent  with  the  rigid  finite  elements  of 
which  it  is  composed.  The  dispersion  relation  for  the  free  surface 
given  by  Eq.  (1)  is  clearly  not  acceptable  in  this  case. 

We  assume  that  the  fluid  slab  is  incompressible  and  inviscid  and 
the  region  z < -d  is  an  extremely  low  density  fluid  so  that  hydrostatic 
pressure  effects  there  are  negligible.  Consistent  with  the  assumption 
of  zero  viscosity,  we  take  the  slab  fluid  motion  to  be  irrotational . 
III.  Analysis 

Since  the  flow  field  is  irrotational,  7 x 3.  = 0 , the  velocity 
field  is  derivable  from  a potential 


£ = 7c(x,z,t) 

Fluid  incompressibility  (2  = constant)  requires  that 

7 . £ = 0 

or  7^  o = 0 


(2) 


(3) 

(^) 


Therefore,  the  velocity  potential  satisfies  Laplace's  equation  within 
the  slab.  Equation  (^)  must  be  supplemented  by  the  associated  boundary 
conditions  in  order  to  solve  for  CD.  At  z = 0,  rigid  wall,  the  vertical 
component  of  velocity  must  vanish 


Acd 

?z 


(z  = 0) 


= 0 


(5) 


At  the  lower  surface,  z = -d,  we  must  proceed  with  care.  The  kinematic 


T 


L 


. . i.:..iarv  la  leza:'::.!;.-; a by  r ;.a;ial,.,;  bcua.aai’y  p>,iata  re- 

l>- 'j:'. -lax’y  a-lr.ca.  Therefore  If  z - -d  is  the  form  of  the 

displaced  lower  surface, 


0 


'.'here  ^ is  the  n.aterial  derivative  defined  by 


it  " d* 


az 


(•') 


seccr.d  crier  terus  in  '^,■3  we  have  from  (c) 


.. 


at  " 


at  z r:  -a 


The  second  boundary  ccnditicn  to  be  specified  at  the  surface  is  the 

* 

dp-T.ar.ic  ccnditicn.  The  finite  surface  eler.ents  separate  the  slab  frcn; 
a constant  cressure  reservoir.  If  an  individual  surface  element  has 
a n.ass,  1::.  then  the  vertical  ncir.ent-uir.  equation  for  each  of  the  eler.ents 
taltes  the  fcLic.ving  forr. 


-f  (p,  - p)  - 


Without  any  loss  of  generality,  we  may  take  the  reservoir  pressure, 
p =0.  This  pressure  merely  establishes  a reference  level  for  the 
entire  pressure  field. 

To  determine  the  pressure  field  fret,  the  velccitp’’  fiela  we  appiy 
the  integrated  form  cf  Zern._ulli^s  ecuatten  tcr  unsteady  irre  ..ataonal 


’If- 


: c 


(10) 


li.ij  ia  valid,  every-ihero  in  the  .:Iab.  In  particular,  applying 
it  at  the  fluid  side  of  the  finite  elements  z = l(x,t)-d  and  neglecting 
second  order  terms,  consistent  with  a linearized  treatment,  we  obtain 


(z  ^ -d)  - f ^ £':(x,t)  = 0 


(i; 


rr.ereiv^re , 


(; 


applied  at  z 

In  o we  Iiave  the  follcwing  equations  to  be  solved  for 

the  dab  ^ z s -d.  The  differential  euuaticn  fcr  the  velocity  po- 
tential iz  Eq.  (h) 


- 0 


(1: 


with  the  bcundiary  ccnditicr.s 


= 0 at  z = 0 
c - 


(1' 


^ n cl 


^9  _ dl(x.t) 
$iz  “ ^t 


(15^ 


at  z = d (15 


The  zclutior.  to  Eq.  (it)  v/hich  satisfies  boundary  condition  (l^+)  a; 


::e:. 


(t)  sin  (k  x) 


B (t)  cos  (k  x)'cosh(k  z) 
n ^ 


V.e  cc;.jiJ.er,  for  the  rxinent,  J^'T'jtecric  perii.dic  initial  atate^: 

- = "(n,0)  periodic  (twice  the  finite  curface  element  length) 


iepicted  in  Fi'.  3«  The  Fourier  nepnecent 


anitiaxicaticn  cc 


•:(x,0)  . 


f:_  i (^) 


3y  c;;,m’j:'.etry,  th.e  full  time  dependent  surface  perturbaricn  muct  remain 
a cpmj;.ctric  cquare  wave  given  by 


i 5 > 


Frcr.  £qs.(l5a)  (ic),  and  (l8)  we  identify 


rrv 

k = , B = 0 

n ^ n 


and  therefore 


= (:<,. ,t;  yt)  (^-i)  =csh 


n=lo,  5 


is  deCermined  from  Eq.  as 


(t'  rr  - - 


\/-\  > 


::.v:uaric 


(20  ai'd  (21,  therefore  represent  a solution  to  Eq.  (ly) 


oatioiiea  tr.e  'OLV.r.aary  CLr.aatacn^  r.0J.  (1-)  and  (l^al.  We 
aiiiticn  .I_y  satis:’;, - the  d;,T;a::;ic  cour.sary  s^r.iitic;: . Es.  (Ipb). 
Fi  iT.l;'.;’  tl'.e  integrand  of  Eq.  (Ipb)  we  have 


¥-i 


£1 


*0  *0  t I'.r* * — r ti.jLG  o6i*i.‘3J  tr^rr'.  bv  t.6!T'*  v;6  r*^'^cv’'^r 

la:,  _ r'  growth  rate  dispersion  relation 


I— -N 


for  the  fc'aitavicr  cf  (t  . This  ir.plies  a grcwth  rate  which  is  de- 


pendent up>.n  th.a  partisular  harrr.cnic  wave  number  each  har- 

monic evolved  with  a growth  rate  given  by  Eq.  (2p),  a square  wave 
initial  profile  would  not  remain  a square  wave. 

Appl.yirwt  the  full  integral  fcimi  cf  dynamic  condition,  Sq.  (Ipb), 
i 'h  talances  the  vertical  force  and  acceleration  for  each,  surface 


^ner»*^  v.**' 


The  growth  rate  scales  as  [V  like  the  Ray le igh-Taylor  rate  corres- 
ponding to  the  finite  element  wave  length,  , The  finite  element  mass, 
m,  is  seen  to  reduce  the  growth  rate  inertially  as  we  would  anticipate. 
Assuming  the  finite  elements  to  be  of  a material  of  density  c * and 
th ickness  d * , 

m = 0 » d * e (27) 


or 


(28) 


For  an  arbitrary  surface  initialization  discretized  by  finite  elements 
of  size  -2,  it  is  anticipated  that  the  maximum  growth  rate  would  scale 
as  demonstrate  this,  below  we  consider  a few  examples. 

For  an  initial  symmetric  surface  perturbation  as  indicated  in 
Fig.  - where  each  wave  length  (X)  corresponds  to  three  pistons  we  would 
anticipate  a maximum  growth  rate  oj  corresponding  to  X/5  = ^ which  would 
scale  as 


(29) 


When  each  wave  length  corresponds  to  - pistons  as  indicated  in  Fig.  5, 
we  obtain  a maximum  growth  rate  of 


~ ® (X/2) 


(?0) 


as  by  symmetry  it  is  not  possible  to  populate  disturbances  of  wave- 
length 2/  . 

The  first  non-trivial  situation  associated  with  the  previous  re- 
sults is  for  the  case  where  the  wavelength  of  the  initial  perturbation 
'•  = '"j?  . This  case  is  illustrated  in  Fig.  6 where,  in  addition,  it  is 
shown  chat  the  surface  displacement  at  any  time  may  be  represented  as 
the  sum  of  two  symmetric  square  waves  of  wave  lengths  2H  and  61, 
respectively. 


10 


That  is,  the  surface  perturbation  is 


(X,  t)  = - 


- ^ nrx 


1 ' 


1 . "StTT'X 

- Sin  ‘ — : — 
n I 


As  before,  the  velocity  potential  in  the  region  0 ^ z ^ -d  satisfies 


7-0  = 0. 


The  boundary  conditions  are 


c;  = o') 

2 /(z  = 0)’ 


t z 


/ (c:.  + gc  )dx  = c) 

J *"  2 . / (z  = — 


Two  dynamic  conditions,  Eq.  (55c)  and  (53d)  are  now  required  to  specify 
the  surface  pressure  balance  for  the  finite  elements.  We  have  con- 
sidered only  the  case  of  massless  pistons  here  (m  = 0).  Proceeding  as 
before  the  solution  for  the  velocity  potential  which  satisfies  the 
differential  equation  (52)  and  boundary  conditions  (55a)  and  (35b)  con- 
sistent with  the  displacement  of  Eq.  (31)  is  given  as 


11 


! for  Che  fundamenCiil  and  u£  ~ f-  ~ for  the  smaller  wavelength  square 

^ Z L 

wave.  Each  above  growth  rate  is  seen  to  scale  in  a Rayleigh-Taylor 

I 

! sense  (Eq.  (1))  with  its  wavelength. 

IV.  Conclusions 

From  the  preceding  analyses,  it  is  seen  that  the  use  of  a plurality 
of  free*pistons  to  drive  a liquid  liner,  while  eliminating  high  fre- 
quency Rayleigh-Taylor  instability  at  the  outer  edge  of  the  liner,  is 
still  subject  to  low-frequency  positional-instability.  The  growth  of 
variations  in  piston  positions  occurs  at  rates  similar  to  those  that 


A 


12 


would  be  calculated  based  on  the  classical  Rayleigh-Taylor  analysis 
for  wavelengths  equal  to  multiples  of  twice  the  piston  size;  signi- 


ficant quantitative  differences,  however,  do  exist.  The  result  of  such 
instability  is  to  provide  a nonuniform,  asymmetric  distribution  of  fluid 
momentum,  which  can  have  serious  consequences  for  both  the  quality  of 
the  inner  surface  implosion  and  the  mechanical  behavior  of  the  rotating 
liner  implosion  system. 
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Fig.  5 — Symmetric  initial  surface  displacement 
for  A = 4t 
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